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\ Abstract: We discuss the synergy of the cosmic shear and CMB lensing experiments to si- 

I multaneously constrain the neutrino mass and dark energy properties. Taking fully account 

CN ■ of the CMB lensing, cosmic shear, CMB anisotropies, and their cross correlation signals, 

' we clarify a role of each signal, and investigate the extent to which the upcoming observa- 

. tions by a high-angular resolution experiment of CMB and deep galaxy imaging survey can 

. tightly constrain the neutrino mass and dark energy equation-of-state parameters. Includ- 

ing the primary CMB information as a prior cosmological information, the Fisher analysis 
reveals that the time varying equation-of-state parameters, given by 'w{a) = wo + Wa{l — a), 
^ ' can be tightly constrained with the accuracies of 5% for wo and 15% for Wa, which are 

■ comparable to or even better than those of the stage-III type surveys neglecting the effect of 

massive neutrinos. In other words, including the neutrino mass in the parameter estimation 
would not drastically alter the Figure-of-Merit estimates of dark energy parameters from 
the weak lensing measurements. For the neutrino mass, a clear signal for total neutrino 
mass with ~ 0.1 eV can be detected with ~ 2-a significance. The robustness and sensitivity 
of these results are checked in detail by allowing the setup of cosmic shear experiment to 
vary as a function of observation time or exposure time, showing that the improvement of 
the constraints very weakly depends on the survey parameters, and the results mentioned 
above are nearly optimal for the dark energy parameters and the neutrino mass. 
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1. Introduction 

Cosmological observations in the last decade have successfully led to the establishment 
of standard cosmological model, and the energy composition in the Universe is well- 
determined with an accuracy less than ~ 10% level. In particular, the observations of 
cosmic microwave background (CMB) and large-scale structure confirmed that our uni- 
verse is described by a flat A cold dark matter (ACDM) model, in which the Universe is 
filled with ~ 30 % of cold dark matter and baryon, and ~ 70 % of the dark energy |5|. 
Together with the first discovery of the late-time cosmic acceleration by the observations of 
Type la supernovae ^] , we are sure that the Universe at the present day is dominated 
by the unknown energy components. 

Taking seriously account of our current understanding and lack of our knowledge of the 
Universe, the cosmology in the coming decade should focus on more profound or even more 
advanced issues; origin and nature of the dark energy and dark matter, and in other words, 
a more accurate description of our Universe beyond a flat ACDM model. The answers to 
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these questions are important, and would shed hght on new physics beyond our current 
knowledge. 

In this paper, we consider the prospects of upcoming observations for tightly constrain- 
ing the nature of the dark energy and the neutrino mass. The theoretical understanding 
of the nature of the dark energy is still limited, and the cosmological observations are the 
only way to reveal the dynamical properties of the dark energy. For this purpose, in the 
literature, the dark energy are often parametrized with its equation-of-state, i.e., p = wp, 
allowing a time dependence of w. Combining several complementary observations, one tries 
to determine and/or constrain tt; as a function of time. On the other hand, determination 
of the neutrino mass is one of the important subject in elementary particle physics, and 
is the key to understand the physics beyond the standard model of particle physics. In 
cosmology, the massive neutrinos play a role of hot dark matter, and they can alter the 
cosmic expansion history as well as the growth of large-scale structure |3^. An inter- 
esting point here is that the cosmological observations can give a tight bound on the sum 
of neutrino masses, and have currently achieved the constraint ^m,/ < 0.2 — l.OeV (95% 



C.L.) (e.g., Refs. |42, 45, ^), which is comparable or even stronger than the constraint 
obtained from the ground-based experiment using /3-decay [^]. Thus, a pursuit of con- 
straining or measuring neutrino mass from upcoming observations is scientifically fruitful 
and at least complementary to ground-based experiments. 

Here, we are specifically concerned with the weak lensing effects on the CMB and 
galaxies, as representative cosmological probes for dark energy properties and neutrino 
mass. The main advantage of the weak lensing effect is to directly map the large-scale 
structure, and we can simultaneously extract the information on the cosmic expansion and 
structure growth imprinted on the lensed CMB map and galaxy images. In particular, the 
weak lensing of galaxies, referred to as the cosmic shear, is expected to be powerful to 
constrain the dark energy and total neutrino mass from future observations, and several 



works on the forecast study have been done |4£, O, 0, H, S^, 13 1. On the other hand, weak 



lensing of CMB, often quoted as CMB lensing, is another sensitive probe of the matter 
density fluctuations at higher redshift, and a signature of lensed CMB has been recently 
detected with a high significance through cross correlation with galaxy clustering |T7|| , 
and directly through the CMB power spectrum ||9|, Future CMB observations such as 
Planck will definitely detect CMB lensing, and a synergy with other cosmological probes 
would help to give a tight constraint on dark energy and neutrino mass. 

As a representative example, we discuss a prospect and a synergy between specific 
upcoming CMB lensing and cosmic shear experiments; cosmic shear survey with the Hyper 
Suprime-Cam (HSC) mounted on Subaru telescope Q, and the CMB lensing experiment 
with the Planck and Atacama Cosmology Telescope with new polarization sensitive receiver 
(ACTPol) ||3^. The cosmic shear survey with HSC plans to observe galaxy images over 
the sky of 2, 000 square degree, with 35 galaxies per square arcminute of the mean redshift 

~ 1. On the other hand, ACTPol has enough sensitivity to measure the CMB lensing 
even at very small scales with higher signal-to-noise ratio than Planck. An important 
characteristic in these two experiments is that survey regions of these observations will be 
overlapped each other, capable of dealing with cross correlation study. 
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0.022 


0.13 


0.72 0.96 2.4 


0.086 0.1 eV 


-1.0 0.0 



Table 1: Fiducial values of cosniological parameters used in this paper. These values are favored 
in the WMAP results l28l. We assume the flat ACDM model with three massive neutrinos. 



Motivated by these situations, in this paper, we address the feasibility and complemen- 
tarity of the upcoming lensing experiments. Based on the Fisher analysis, we quantitatively 
investigate how well we can constrain the dark energy and total neutrino mass, taking fully 
account of the CMB lensing, cosmic shear, CMB anisotropics, and their cross correlations. 
To our knowledge, this is the first work to consider full cross correlation signals between 
the two lensing signals, including the effects of both the dark energy and the neutrino mass. 
There have been several forecast studies on both the CMB lensing and cosmic shear, but 
these are restricted to specific models such as the early dark energy |l^ or those neglecting 
the massive neutrinos [ 20 1 . We here discuss a role of each signal and the influence of param- 



eter degeneracies in constraining the neutrino mass and dark energy properties. Further, 
dependence of the forecast results on the survey setup is investigated in detail. In this 
respect, the analysis in this paper includes a general result not only valid for the specific 
lensing surveys with HSC and ACTPol, but also applicable to other lensing experiments. 

This paper is organized as follows. In section |2[ we briefly review the CMB lensing and 
cosmic shear. In section |^, we summarize the Fisher matrix formalism, and describe the 
canonical setup of the CMB lensing and cosmic shear experiments with Planck, ACTPol 
and HSC. Then, in section ^, the signal-to-noise ratio and forecast results for parameter 
estimations are presented in detail, focusing on the dark energy and total neutrino mass. 
In section ||, to see the sensitivity of the forecast results to the details of the survey setup, 
we allow to vary the setup of cosmic shear experiment, and examine how the constraints 
on neutrino mass and dark energy properties are changed. Finally, section ^ is devoted to 
summary and conclusion. 

Throughout the paper, we calculate the power spectra for a fiducial set of cosmological 
parameters in Table i.e., the density parameter of baryon i^^^h"^, of matter Q^h"^, dark 
energy density fl^, scalar spectral index rig, and scalar amplitude Ag at k = 0.002 Mpc^^, 
reionization optical depth r, total neutrino mass ^ mj^, and dark energy equation-of-state 
parameters, wq and Wa, with the functional form of w{a) = wq + Wa{l — a) (e.g., pO[|). 

2. Weak gravitational lensing of CMB and galaxies 
2.1 Power spectra of vi^eak gravitational lensing 

Primary CMB temperature, 0(n), and polarization fields, E{n) and B{n), are distorted 



due to the large scale structure between the last scattering surface and us (see |32| for 
a review). The weak lensing effect on the CMB fields is characterized by the deflection 
vector, d such as 

e(n) = e(n + d), (2.1) 
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where Q denotes the lensed temperature field, and the same is as well with the case of 
polarization fields. Note that lensed i?-mode polarization can be generated by lensed 
distortion of primary E-mode |58|. The deflection angle due to the gravitational lensing is 
obtained by solving the geodesic equation, yielding 

f'X* 

d{n) = -2 dx^^^^Vf,ij{xn,rjo-x), (2.2) 
Jo X*X 

where x* denotes the comoving radial distance to the last scattering surface, and r/o is the 
conformal time at present, and ijj{r,r]) is the Newton potential. Note that we set units for 
the light velocity. The deflection field thus traces the gravitational field produced by the 
structure between the last scattering surface and us. Using Eq. (|2.2|) , a statistical quantity 
of our interest, i.e., the angular power spectrum of the deflection angle, Cf'^, is theoretically 
computed as [20| 

Cf = l I ^nMk)[A^{k)]^ (2.3) 
where Pinit(^) is the matter power spectrum at an early time. The function A"'(fc) is given 

by 

V (^-1)! Jo X*X a[x) 

with $7m and -f^o being the density parameter of mass and the Hubble parameter at present 
time, respectively. The function D{k] z) represents the growth factor, defined by the square 
root of the ratio of matter power spectra, D{k;z) = \J P{k; z) / Pinitik) ■ Note that the 
scale-dependent growth naturally arises from the non-linear gravitational evolution and 



free-streaming suppression by the massive neutrinos (see 2.2) 



The deflection fleld can be efficiently reconstructed from the non-Gaussian nature in 
the lensed CMB statistics |21, In the flat-sky approximation, the estimator for the 
deflection field is given by 

d^yH) = I -0^X(i')Yie-e')gxYii,£'), (2.5) 

where the subscripts X and Y denote the lensed temperature and polarization fields, i.e., 
B, E, B, and {£) and gxvif-^f-') denote the normalization and the optimal weight 
function, respectively (see for explicit expressions for A^^ {I) and gxvit-, (■'))■ The 
essence of this quadratic estimator is that information from different scales are used to 
reconstruct the defiection field. The optimal weight is chosen so as to minimize the lensing 
reconstruction error, Nf'^. If one consider the temperature field only, for an illustrative 
example, the error Nf'^ is estimated as 



(2^)2 2(Cfe + Aree)(cee , ^ee ) 



(2.6) 



where the quantity Nf® is the temperature noise power spectrum (see section ^). In this 
paper, we compute Nf'^ following the technique developed in ||3^ in which the temperature 
and polarization fields are optimally combined in the full-sky treatment. 



- 4 - 



On the other hand, in a weak lensing survey of galaxies, we first measure the elhpticity 



of each galaxy image, and then estimate the shear field 7(ri) [36|. Although the shear field 
is projected onto two-dimensional sky, the inclusion of redshift information significantly 
improve the sensitivity to the cosmological distance and structure growth, which leads to 
a tighter constraints on the dark energy. This technique is referred to as the tomography 
I IS], and with the photometric redshift samples divided into several subsamples binned 



with redshifts, we obtain various combination of cross shear spectrum, e.g., Cj^'^^ for cross 
spectrum between z-th and j-th redshift bins. Given redshift distribution of source galaxies 



in i-th, ni{z), the shear power spectrum is calculated from ||20[] 

=IJ ^Pinit{k)A^^{k)A^^{k), (2.7) 



where Aj^{k) is given by 



A..(A-) = ,m±^^ r,.,mr'<IX^^^^^¥^i.m. (2.8) 



The quantity nj is the average number density per square arcminute in i-th bin, defined as 

;>oo 

fii = I ni{z)dz. (2-9) 



JO 

Note that apart from the prefactor {i+2){i—l)/A, the expression ( p.7D with ( |2.8[) is reduced 
to that of the CMB lensing ( |2.3| ) if we choose n{z) = 5d{z — z^,), where (J^j represents the 
Dirac's delta function and means the redshift at the last scattering surface. Thus, the 
shear power spectrum Cj'^ generically has a larger power at small scales {(. ^ 1) than the 
CMB lensing spectrum Cf'^. 

In addition, for a survey region overlapping between the cosmic shear and CMB lensing 
experiments, the cross correlation signal between the deflection and shear fields can be 
obtained. Since both the CMB lensing and cosmic shear signals are originated from the 
large-scale structure, the cosmological information from the two lensing measurements is 
statistically correlated. Therefore, the cross correlation signal must be properly taken into 
account in cosmological analysis. The cross power spectrum between the CMB deflection 
and galaxy shear flled is predicted as [^] 

Cp' = 11 ^Pinit{k)A^{k)A''^{k). (2.10) 

Finally, another important cosmological signal is obtained through the late-time in- 
tegrated Sachs- Wolfe (ISW) effect, which results from the late-time decay of the gravita- 
tional potential fluctuations induced by the departure of cosmic expansion from that of 
the Einstein-de Sitter universe. The ISW effect produces a non-vanishing cross correlation 
between the CMB temperature and weak lensing fields, and its characteristic signal ap- 
pears at large angular scales. The angular power spectrum of CMB temperature and weak 
lensing field X (i.e., X = d or -ji) becomes ||2^ 



2 f dk 

= -J ^^'init(A:)Aisw(fc)A^(fc), (2.11) 



- 5 - 



where the kernel Aisw(^) is given by 

Aisw(A;) = 3nmHl£ dz^ {^^^} •?'^(^^(^))- (2.12) 

All the spectra presented above include information on the evolution of matter power 
spectrum in a somewhat complicated manner. In next subsection, we will discuss the cos- 
mological parameter dependence of these power spectra, especially focusing on neutrino 
mass and dark energy properties. Throughout this paper, all the power spectra are com- 



puted with the modified version of CMB Boltzmann code, CAMB |33], with the fiducial 
cosmological parameters given in Table ||. For nonlinear matter power spectrum relevant 
for the auto and cross power spectra of cosmic shear and deflection fields, we adopt the 
fitting formula proposed by Smith et al (2003) [^] . 

2.2 Cosmological information from CMB lensing and cosmic shear 

To see how the lensing spectra depend on the cosmological parameters, in Figure |l], we plot 
the logarithmic derivatives of angular spectra Cf'^, C'J^'^^ and C^'^^ with respect to the dark 
energy equation-of-state parameters wq, Wa (left), and the neutrino mass characterized by 
the dimensionless parameter, = (^ mi,/93.14eV)/r2m/i^ (right). Here, the quantity 73 
represents the shear field estimated from the galaxy subsample in the third redshift bin, 
whose precise meaning will be specified later (see Sec. |^. Note that the cosmological 
dependence of the other shear fields from different redshift bins is qualitatively the same. 
For comparison, we also show the dependence on matter density, i.e., the logarithmic 
derivative dlnCi/dlnQ^h'^ . 

The influence of dark energy on the lensing power spectra basically comes from the 
growth factor of matter fluctuations, D{k;z), and because of the homogeneous nature of 
the dark energy, its effect on growth factor is nearly independent of scales. Thus, the 
resultant variation of angular spectra with respect to the dark energy equation-of-state 
parameters is almost scale-independent. On the other hand, massive neutrinos lead to a 
characteristic suppression on the growth of large scale structure below the free-streaming 
scale. Since the characteristic wavenumber of the free-streaming scale becomes large as 
increasing the neutrino mass, as a consequence, the logarithmic derivative of lensing spectra 
can be negative, and decreases as increasing multipoles. That is, the amplitude of lensing 
spectra is reduced with the effect of neutrino mass at small scales, and with a high angular 
resolution lensing experiment, we can detect a clear signature of free-streaming suppression. 

These scale-dependent or scale-independent natures of the lensing power spectra are, in 
principle, very powerful to constrain the neutrino mass and dark energy. However, we note 
here that there exist some parameters that exhibit a similar scale-(in)dependence, which 
can be the source of parameter degeneracy. As shown in Figure |^, the logarithmic derivative 
of the lensing spectra with respect to the quantity In fim/i^ gives a similar trend to that of 
the neutrino mass parameter, dlnC^/ dfy, and thus the density parameter Q.^h? can mimic 
a scale-dependent suppression by massive neutrino. Also, the equation-of-state parameters 
wq and Wa can put a similar scale-independent behavior in the lensing spectra, and it seems 
difficult to discriminate between each other even if we use the lensing information at small 



-6- 



-1.2 '^^3 

- - - 7373 

10 100 1000 

I 

Figure 1: Logarithmic derivatives of the lensing power spectra, Cf"^, CJ^''^ and Cf^^ ^ with respect 
to Wo and Wa (left), and fu = mi,/93.14eV)/f2in^^ and ^^h? (right). The derivatives are 
evaluated around the fiducial set of cosmological parameters summarized in Table |l|. To compute 
the derivative of cosmic shear spectra, the galaxies are divided into the three redshift subsamples 
as described in Sec.^, and we especially plot the results for the third subsample (i.e., 1.5 < z). Note 
that to see similarity of the parameter dependence, the result of c? In C^/c? In fim/i^ is multiplied by 
-0.45. 
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scales. Note that the primary CMB information alone is insensitive to the dark energy 
properties and it exhibits some degeneracies, e.g., between InOm/i^ and f^A. Hence, even 
combining the CMB data, there may remain a sizable amount of degeneracies among wq and 
Wa, 'Yl'^f ^'^d Q^h?. In this sense, the degree of improvement for parameter constraints 
seems rather non-trivial. For more quantitative aspect of the parameter estimation, we 
will proceed to the Fisher analysis. 



3. Fisher matrix formalism 



In this section, we summarize the Fisher matrix formalism used in the subsequent analysis, 
and describe the canonical setup for CMB and cosmic shear experiments, for which forecast 
results are presented in section ^ 

Given the angular power spectra theoretically parametrized by a set of parameters p, 
the cosmological information on these parameters obtained from the combination of several 
experiments is quantified by the Fisher matrix [55|: 



i^.- = E ^/skyTr , (3.1) 

where the quantity Ci represents the covariance matrix for angular power spectra given 
below. Pi is a cosmological parameter which we want to estimate, and the /sky is the sky 
coverage. Using the Fisher matrix, the 1-a (68%C.L.) statistical uncertainties for the cos- 
mological parameter pi marginalized over other parameters, cr{pi), is estimated as cr{pi) = 
\J {i^^ijjj. Also, the statistical correlation or degeneracy between parameters pi and pj 
can be deduced from the off-diagonal component of inverse Fisher matrix and is 

quantified by defining the correlation coefficient, r{pi,pj) = {F~^}ij /{{F~^}ii{F~^}jj}^/'^ . 



It is worth noting that Eq.( |3.lD rehes on the assumption of the Gaussianity for hkehhood 
function, which is incorrect in some situations. In particular, cosmic shear is sensitive to the 
nonlinear gravitational evolution, which leads to the non-Gaussian covariance. However, 
the actual impact of non-Gaussian covariance on the parameter estimation is turned out to 
be insignificant, and it degrades at most few percent level [^, 52|. Although this does not 
imply the validity of the error estimation with Gaussian likelihood function (for example, 
in the case of neutrino mass, the condition, ^rriu > 0, leads to non-Gaussian error pC|]), 
we adopt Eq.( |3.lD to investigate the potential power of weak lensing experiments. 

Throughout the paper, we specifically consider Planck, ACTPol and HSC as the repre- 
sentative CMB and cosmic shear experiments, from which we can obtain the temperature 
(0) and (E-mode) polarization {E) data for the primary CMB anisotropies, the deflection 
angle (d) data for the CMB lensing, and the shear field data of galaxies divided into several 
redshift bins (7^; i = 1, • • • ,N). Then, the full covariance matrix, Ci, is written in the 
form as 
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(3.2) 



(^7n7n ^^7n7ny 

Here, is the noise power spectrum. The amplitude and shape of the noise spectra 

depends on the survey design for the CMB and lensing experiments, which will be 
discussed below. 

For the CMB experiment, Planck, we use seven frequency channels for temperature 
and polarization observations, assuming the sky coverage of /sky = 0.65. As for the ground- 
based experiment, ACTPol, single frequency channel with = 148GHz is used for high- 
resolution observation, and we assume /sky = 0.1 to enhance the CMB lensing information^. 
When estimating the temperature and polarization power spectra from ACTPol, we com- 
bine the Planck data of i < 700, in order to remedy a large uncertainty at large angular 
scales arising from the atmospheric temperature fluctuations. 

In both Planck and ACTPol, the dominant noise source for temperature and polar- 
ization observations may be the photon shot noise. The noise power spectra are then 
expressed as 



XX 



XX\-1 



N, 



XX 



\TcMB J 



81n2 



(3.3) 



with TcMB = 2.7K being mean temperature of CMB. Here, the quantity 6^^ is the beam size, 
and cTjy represents the sensitivity of each channel to the temperature cTu^t or polarizations 



^This setup is somewhat different from the original proposals of Ref. jS^, but we here keep it to inves- 
tigate the potential power of high resolution experiment. 
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Experiment /sky v [GHz] Qy [arcmin] Gy^T [/xK/pixel] a^^p [/iK/pixel] 



Planck @ 065 30 33 44 

44 23 6.5 9.2 

70 14 9.8 13.9 

100 9.5 6.8 10.9 

143 7.1 6.0 11.4 

217 5.0 13.1 26.7 

353 5i} 4(U 81.2 

ACTPol 11 Ol 148 L4 3^6 slT 



Table 2: Assumed experimental specifications for the Planck and ACTPol. The quantity is the 
beam size, and Oi, represents the sensitivity of each channel to the temperature g^^t or polarizations 
CTiy.p, depending on the power spectrum of temperature (X ~ 0) or polarizations (X = E or B). 
The quantity v means a channel frequency. 



(Tj/^P, depending on the power spectrum of temperature {X = Q) or polarizations {X = E 
or B). Specific values for these quantities are summarized in Table Note that the 
foregrounds contamination in the CMB data would be additional noise source for lensing 
experiments, and it could not only degrade but also bias the cosmological constraints. 
Recent estimate by Ref . suggests that the foregrounds contribution to the angular power 
spectrum would become less than the instrumental noise, if foregrounds can be successfully 
subtracted at a 1% level. Although this is still challenging and a more elaborative study is 
necessary for the foreground removal, we here ignore the effect of foreground contamination 
in order to explore the potential and complementarity of the lensing experiment. 

For the CMB lensing data, photon shot noise given by Eq. ( p. 3D is also the dominant 
noise contribution, but this time, we must further consider the errors associated with 
reconstruction technique. We adopt the optimal quadratic estimator proposed by Okamoto 
&: Hu (2003) [^, and the temperature, and E- and B-mode polarization data are used to 
estimate the power spectrum of deflection angle. The noise power spectrum N^'^ is then 



computed based on the expression (42) in Ref. |38|. 

On the other hand, canonical setup for cosmic shear survey with HSC discussed here 
roughly match the survey plan proposed in Ref. Q. We consider the deep imaging survey 
with area 2,000 deg2 and mean redshift = 1, assuming the redshift distribution of 
galaxies given by 

3N_„ 
2z^ 



n(z) = — |- exp 



"0 



1.5 



(3.4) 



with the mean number density of galaxies, Ng = 35 arcmin"^. Note that zq is related to 
the mean redshift through zq = 0.692^. We divide the whole galaxy sample into the three 
redshift subsamples (i.e., A'^ = 3) for lensing tomography, and use the auto and cross power 
spectra between different redshift bins; < z < 0.7, 0.7 < z < 1.5, and 1.5 < z. Note 
that for lensing tomography based on the photometric redshift technique, the uncertainty 
arising from the photometric redshift error is crucial for the cosmological analysis. To 
mimic this effect, we suppose that the photometric redshift estimates are distributed as a 
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Table 3: Specific values of Api used to evaluate the derivatives of angular power spectra, Eq.(3.10) 



Gaussian with rms fluctuation cr{z). Then the actual redshift distribution for i-th galaxy 
subsample becomes p^] 



r I Zi — Z \ I Zi^l — Z 

erfc —= — erfc 



(3.5) 



,V2a{z)J \V2a{z) 
where erfc(x) is the complementary error function defined by 

erfc(x) = / dz exp{—z'^). (3-6) 

Jx 

For simplicity, we adopt the scaling relation for the photo-z error: 

o-(z) = 0.03(1 + (3.7) 

Apart from the calibration systematics for shear estimation, the main noise source for 
cosmic shear measurement is the intrinsic ellipticity of galaxies, which can be described as 

Nr'=6ij^4^. (3.8) 

The quantities (7int) and hi are the rms intrinsic shear and the number density of 
galaxies per steradians in the i-th bin, respectively. We adopt the empirically derived 

1 /2 

value, (7int) = 0-4 0- Note that the quantity hi is related to the average number 
density per arcminute square in the i-th bin, n^: 

180 \^ 



hi = 3600 nj ( — J str-^ (3.9) 



Finally, in computing the Fisher matrix, Eq.( |3.lD , we replace the derivatives of the 
angular power spectrum Ci with respect to the cosmological parameters pi with a finite 
difference given by 

_ CejPi + ^Pi)- Cijpi - Api) 
dpi 2Api 
Specific values of the difference Api used in the analysis are listed in Table ^. 

4. Results 

In this section, forecast results for power spectrum measurements and parameter estima- 



tions are presented based on the canonical setup in previous section. In section 4.1, we first 
compute the signal-to-noise ratio for CMB lensing and cosmic shear observations. We then 
present the results of Fisher analysis for parameter forecast in section just focusing on 
total neutrino mass (Ylifny) and dark energy equation-of-state parameters {wf} and Wq)- In 



section a role of cross correlations on the parameter constraints is discussed in some 
details. 
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4.1 Signal-to- noise ratio 

To see how robustly the lensing power spectra can be measured with a high significance 
at each scale, let us first look at the expected signals and noises from the CMB lensing 
and cosmic shear measurements. Figure I plots the power spectra for CMB lensing (top 
left), cosmic shear (bottom left), and their cross correlation (bottom right) for fiducial 
cosmological model in Table Also, we plot the variants of power spectra with slightly 
different values for wq and ^m^. The plotted errors for Planck, ACTPol, and HSC are 
estimated from 



AC, 



XY 



V(£+l/2)/skyA^ 



{21 + l)/skyA£ 



{X = Y) 



(4.1) 



where /sky is the sky coverage of each experiment, and A£ is the size of multipole bin, for 
which we set A£ = 200. In Figure ^, we also show the signal-to-noise ratio for each power 
spectrum given by 



{S/N) 




r<xY 



XY 



(4.2) 



Note here that we set = 1 to evaluate AC^-^. 

In general, cosmic shear signal becomes larger and has a higher signal-to-noise ratio at 
smaller angular scales, and for the higher redshift sources. Since the free-streaming sup- 
pression of the massive neutrinos is known to appear at relatively smaller scales (or higher 
multipoles), the cosmic shear signals are generally sensitive to the change of the total mass 
of neutrinos, as well as to dark energy properties through the late-time variation of struc- 
ture growth. Figure |2| indicates that the cosmic shear signals from HSC survey potentially 
have enough sensitivity to detect total neutrino mass of ~ O.leV and to constrain dark 
energy properties. 

By contrast, the angular power spectrum of deflection angles from CMB lensing differs 
from that of the cosmic shear, and has a larger amplitude at lower multipoles, as mentioned 
in section ^. Thus, the CMB lensing seems less sensitive to the suppression effect of 
massive neutrinos, and a high-angular resolution experiment is required for measuring 
the neutrino masses with a sub-eV. As shown in Figure |2[ with ACTPol, we can clearly 
discriminate between the total neutrino masses of the difference ~ 0.2eV from the CMB 
lensing experiment alone. Note that the signal-to-noise ratio for CMB lensing from the 
ACTPol is rather higher than that of the cosmic shear signal obtained from the HSC 
survey. We thus naively expect that the ACTPol can give a better constraint on neutrino 
mass than the HSC. In practice, however, there are a sizable amount of degeneracies among 
several cosmological parameters. In particular, as indicated by Figure |^, the neutrino mass 
is tightly correlated with ^^Y? . Thus, the final outcome of the neutrino mass constraint. 
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after marginalizing over the other parameters, cannot be straightforwardly understood only 
from the signal-to-noise ratios. 

Finally, we note that the cross correlation signal d'ji obtained from ACTPol and 
HSC has a high signal-to-noise ratio (bottom-right panel of Figure |2|). We also compute 
the signal-to- noise ratios for other cross correlation signals, and find that temperature- 
deflection cross correlation Qd from ACTPol and Planck would be measured with high 
signal-to- noise ratio, S/N~ 0(10), while the temperature-shear cross correlations Q'ji have 
a low signal-to-noise ratio, S/N< 1. These results basically come from the facts that the 
non-vanishing cross correlations can be attributed to the ISW effect, and a large signal of 
ISW effect only appear at large-angular scales. 

4.2 Parameter forecast 

We now present the forecast results for cosmological parameters derived from Planck, 
ACTPol and HSC. Table ^ summarizes the marginalized 1-a (68%) errors for the results 
combining two (Planck and HSC) or three observations (Planck, ACTPol and HSC). In 
each case, we examine following three cases: 

• +d : including CMB lensing data {Cf'^ and Cf'^) 

• +7 : including cosmic shear data [C]'^' and Cf^') 

• +d + 7 : combining aU power spectra, i.e., Cf'^, Cf'^, C]''^\ Cf"'\ and Cp^ 

Note that we use the lensing information up to £max = 3000, and add the primary (unlensed) 
CMB data (i.e., C®®, Cf^ and Cf^) as prior information in all cases. We assume that 
the observed area of HSC is entirely overlapped with that of ACTPol, and similarly the 
survey region of ACTPol is totally included in the nearly full-sky survey with Planck. 

Overall, the constraints obtained from the single experiments {+d or +7) are almost 
comparable, and combining all the lensing observations {+d + 7) moderately improves the 
constraints on cosmological parameters. The constraints on wq and Wa obtained from the 
Planck and HSC roughly correspond to those from the Stage-HI class experiments defined 
by Dark Energy Task Force Q (e.g., the Dark Energy Survey and Pan-STARRS-4), and 
the constraints are further improved by a factor of ~ 2 when we add the ACTPol data. 
The neutrino mass constraint is also improved by a factor of ~ 1.5, and with Planck, 
ACTPol and HSC, the total mass of neutrinos can be detected with ~ 2-0" significance 
for a fiducial value = O.leV. For other cosmological parameters, the constraints 

by adding ACTPol data are 1.5 — 2 times better than those obtained from Planck and 
HSC. Note that our neutrino mass constraint from the CMB lensing is roughly consistent 
with those obtained by Refs. |l^, 26]. Also, the constraints on wq and Wa from 



cosmic shear measurement roughly match the results obtained in Ref. |T^. In other words, 
including the neutrino mass in the parameter estimation would not drastically alter the 
constraints on dark energy equation-of-state parameters from the lensing measurements. 

In Figures ^ and ^, to elucidate the impact of upcoming lensing experiments on neutrino 
mass and dark energy properties, we plot the expected 1-a (68%C.L.) contours on hiO^^h? 
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Figure 2: Angular power spectra for the deflection angle C^'^ (top left), cosmic shear C/'^^ (bottom 
left) and deflection-shear cross correlations C^^' (bottom right). The expected error bars are 
estimated from the canonical survey parameters with A£ = 200. For cosmic shear survey, the 
observed galaxies are divided into three subsamples; < z < 0.7, 0.7 < z < 1.5 and 1.5 < z. Also 



we show the signal-to-noise ratio {S/N)^£ defined by Eq.(1.2) in the bottom of each panel 



and ^m,y, and wq and Wa, respectively. Left panels are the results obtained from the 
Planck and HSC survey, while the right panels are obtained by further adding the ACTPol 
data. Comparing with the constraints coming from the primary CMB data alone, the 
size of error ellipses becomes rather reduced when we add the lensing information, but 
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Planck + HSC 


Planck + ACTPol + HSC 


parameter 


+d 


+7 


+(i + 7 


+d 


+7 


+(i + 7 


ln(Ob/i^) 


0.0063 


0.0058 


0.0057 


0.0030 


0.0030 


0.0029 


ln(J)m/i2) 


0.016 


0.011 


0.010 


0.0094 


0.0080 


0.0074 


Qa 


0.031 


0.022 


0.021 


0.017 


0.014 


0.013 


Wq 


0.14 


0.12 


0.11 


0.056 


0.054 


0.052 


Wa 


0.28 


0.24 


0.21 


0.18 


0.17 


0.14 


ris 


0.0033 


0.0029 


0.0028 


0.0023 


0.0023 


0.0022 


HA) 


0.014 


0.012 


0.012 


0.011 


0.011 


0.010 


T 


0.0040 


0.0040 


0.0040 


0.0040 


0.0040 


0.0040 


Y^mu [eV] 


0.13 


0.060 


0.055 


0.061 


0.046 


0.041 



Table 4: Forecast results for marginalized 1-cr errors for each cosmological parameter, assuming 
Planck and HSC (left three columns) and Planck, HSC and ACTPol (right three columns). The 
labels, +d, +7, and +d + 7, respectively indicate the cases including CMB lensing, cosmic shear, 
and combining both. The primary CMB information is included in all cases. 

the relative differences between the three cases (i.e., +d, +7 and +d + 7) are basically 
small, as we mentioned above. Further, the inclinations of the ellipses are almost the same. 
This can be partly explained by the fact that cosmological information on dark energy 
and massive neutrinos mostly comes from lower redshifts {z < 1), and thus, the cosmic 
shear and CMB lensing power spectra have almost identical information on the dark energy 
properties and neutrino mass, though the CMB lensing signal has potentially sensitive to 
higher redshifts. As shown in Figure there exists non- vanishing correlation d'yi, and 
thereby the cosmic shear and CMB lensing signals cannot be regarded as statistically 
independent signals. As a result, only a moderate improvement of the constraints on the 
parameters, wq and Wa, is achieved, and a tight correlation between cosmological 

parameters, ^rriy and InVL^h?, wq and Wa, still remains in the combined results of two 
lensing experiments. Nevertheless, the CMB lensing and cosmic shear experiments greatly 
improve the constraints from the primary CMB information, and these can be used as 
an independent cross check for extracting cosmological information in an unbiased way. 
Hence, science benefit for combining two lensing experiments is still valuable. 

4.3 Role of cross correlation statistics 

Here, we discuss the role of the cross correlation statistics to the parameter constraints 
given in Table |^. Apart from the CMB temperature-polarization cross correlation &E and 
tomographic lensing correlations 7i7j, there are three kinds of measurable cross correla- 
tions, i.e., Qd, Qji, and d'ji. Among these, the deflection-shear cross correlations have 
almost identical information to the deflection or shear auto correlation, and they do not 
significantly contribute to the neutrino and dark energy equation-of-state parameters, as 
we already discussed in previous section. Also, the temperature-shear cross correlation has 
a small signal-to-noise ratio, and no valuable cosmological information can be obtained. 
On the other hand, the signal-to-noise ratio for temperature-deflection cross correla- 
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Figure 3: Expected 1-a error contours on Infl^h^-J^'n^u plane assuming Planck; and HSC (left), 
and Planck, ACTPol and HSC (right). The gray dotted lines represent the results ft'om primary 
CMB, while the red solid and green dashed lines indicate the results further including CMB lensing 
and cosmic shear, respectively. The blue filled ellipse are the results of all measurements. 
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Figure 4: Same as Figure |3| but for expected 1-a error contours on wg-Wa plane. 



tion is S/N~ 0(10) for the Planck and ACTPol experiments. Since the non- vanishing 
contribution of this cross correlation basically comes from the ISW effect, it is expected 
to be sensitive to the early-time evolution of dark energy. Although this contribution has 



been previously ignored in the Fisher analysis in Ref. |14], where they compute the con- 
straints on wq and Wa from the CMB lensing signal alone, we find that this signal gives 
an interesting contribution to the constraint on the time dependence of the dark energy 
equation-of-state parameter, Wa- 

To see the significance of the contribution from temperature-deflection cross corre- 
lation, in Figure |5|, we plot the two-dimensional contour of 1-a error on Wa-^A plane, 
assuming the CMB lensing data from Planck and ACTPol. Here, in addition to the con- 
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straints coming from the primary CMB information (dashed, labeled as CMB), we consider 
the following three cases: 

• +dd : using Cf'^, 

• +Qd : using Cf'^, 

• +ed + dd : using Cf'^, Cf 

In all three cases, we add the primary CMB information (C®®, Cf^, and Cf^), and 
in computing the Fisher matrix for +dd and +Qd cases, the contributions coming from 
the CMB lensing auto/cross correlation are separately computed and are then added to 
the Fisher matrix for primary CMB information. Note here that the +Qd + dd case just 
coincides with the +d case examined in previous section. 

As we see from Figure |5|, the temperature-deflection cross correlation can give a slightly 
better constraint than that from the deflection auto correlation. This is somewhat remark- 
able in the sense that the signal-to-noise ratio for temperature-deflection cross correlation 
is rather lower than that of the deflection auto correlation. The reason for this is pre- 
sumably explained by the fact that the inclusion of the cross correlation breaks the de- 
generacy between i}\ and Wa- the correlation coefficient between and Wa, given by 
r(r2A, Wa) = {F~^}ij/[{F~^}ii{F~^}jj]~^/'^ with i = Wa and j = Oa, changes the value 
from r{Q\, Wa) = 0.3 to —0.1, relatively reducing the statistical correlation between these 
parameters. The result indicates that the temperature-deflection cross correlation has a 
better sensitivity to higher redshifts compared to the auto correlation, and potentially 
gives a tight constraint on the time variation of dark energy equation-of-state parameter, 
w{z). In appendix this point will be studied in detail based on the principal component 
analysis. 



5. Sensitivity to the survey design of cosmic shear experiment 

In this section, to elucidate the robustness or sensitivity of the forecast results in previous 
section, we allow to vary the setup of cosmic shear experiment, while we keep adopting 
the Subaru HSC facility. Then, we examine how the constraints on the neutrino mass and 



dark energy properties are changed. In section 5.1, we first discuss how to characterize 
several key parameters of the cosmic shear experiment as a function of observation time or 
exposure time. In section ^.2| , the forecast results for the constraints on neutrino mass and 
dark energy properties are given, and the sensitivity to the observation time or exposure 
time is investigated. 

5.1 Modeling cosmic shear surveys 

In a cosmic shear survey, the exposure time, total observation time, and number of redshift 
bins are essential key parameters to determine a survey design, and these are related to 
the sky coverage, survey depth and number of galaxy samples. Thus, the final outcome 
of the parameter constraints can be sensitively affected by those parameters. Here, we 
wish to relate the parameters of survey design, and try to characterize the forecast results 
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Figure 5: Expected 1-cr error contours on ^t^-Wa plane. The dot-short-dashed hne (magenta) 
show the constraints obtained from deflection auto correlation and the dot- long-dashed line 

(cyan) is the result from temperature-deflection cross correlation (-|-0d). The red solid line is the 
combined results from all signals of the CMB lensing. Note that primary CMB information is 
added in all three cases. As a reference, we also plot the result from primary CMB information 
(gray dotted). 



for neutrino mass and dark energy equation-of-state parameters as functions of the most 
important parameters, i.e., total observation time and exposure time. To do this, we 
basically follow the treatment by Yamamoto et al. |57|. That is, for a given number of 
redshift bin N ^ we relate the sky coverage /skyj mean redshift z^, the total number density 
of galaxies per square arcminute Ng introduced in section ^ to the exposure time iexp 
and total observation time Tobs- Supposing the Subaru HSC facility for a cosmic shear 
experiment, we then express the forecast constraints as functions of Tobs and texp- 

In the tomographic lensing technique, several filters are used to construct the redshift 
subsamples, and the available number of redshift bins is basically determined by the number 
of filters. Thus, we must first specify the number of filter and determine the exposure time 
for each filter. Following Yamamoto et al. |57], we assume that the i-band filter is used 
for the case = 1, and the i- and r-band filters are used for N = 2. As for the cases with 
= 3 and 4, we consider g-^ r-, i- and z-band observations. Based on the survey proposal 
in Ref. for the one field-of-view observation, the exposure time of the r- and g-hands 
is assumed to be equal to that of the i-band filter, but the z-band exposure time is set to 
be 4/3 times longer. Then, depending on the number of filters or redshift bins, the total 
exposure time ^ texp per one field-of-view is effectively expressed as the i-band exposure 
time, which we denote by texp- The relation between ^iexp and fgxp is summarized in 



-17- 



Table |. 

We now relate the parameters Zm, Ng and /sky to the exposure time texp and total 
observation time Tobs- The mean redshift Zm for the observed galaxy distribution depends 
on the survey depth, and is related to the exposure time per one field-of-view. According 
to Ref. [^], a simple scaling relation between z^ and texp is given by 

\ 0.067 

•2m(iexp) = Zrafl I — 7 1 , (5-1) 

V *0 / 

where we set to = 15min and Zmfi = 1^. As for the total number of galaxies per square 
arcminute Ng, we adopt the empirical relation determined by the pilot observation |l|: 

ft 

iVp(texp) = 35 X i-j^j arcmin-2. (5.2) 

Finally, adopting the 1.8 deg^ of the field-of-view of HSC [||], the sky coverage /sky is 
expressed as the function of texp and total observation time Tobs- According to Ref. p7|| , 
we have 

/sky(texp) = 1.8deg2 . X — (5.3) 

i • i lexp + love 2 , 000 deg 

Here, we take into account the overhead time tove = 5 min and the processing time 0.1 x 

^exp • 

To sum up, given the exposure time texp and total observation time Tobs, the redshift 
distribution of galaxies characterized by Eq. ( |3.4[ ) is specified and the sky coverage of cosmic 



shear survey is fixed through the relations (5.1)-(5.3), depending on the number of filters 



or redshift bins (see Table ^). Then, we can proceed to the Fisher analysis just following 
the procedure in section |3|. Finally, we note that the canonical setup for the cosmic shear 
experiment examined in section Q corresponds to the description mentioned above with the 
specific values of texp = 15 min, Tobs = 180 nights^ and = 3. In Figure ^, as a reference, 
the dependence of the sky coverage /sky on the number of redshift bins N is plotted against 
the i-band exposure time texp, fixing the observation time to Tobs = 180 nights. 

5.2 Results 

Based on the procedure in the previous subsection, we vary the setup of cosmic shear 
survey, and derive the constraints on the neutrino mass and dark energy equation-of- 
state parameters, characterized by the exposure time and total observation time. In what 
follows, combining the CMB lensing and primary CMB information obtained from Planck 



and ACTPol, the results are separately presented for neutrino mass (Sec. 5.2.1 ) and dark 



energy properties (Sec. |5.2.2D . 



^Ref. adopted to ~ 30 min and Zm,o = 0.9 to just follow the parameters described in j^. Here, we 
adopt different values to match the recent study on the HSC survey plan [^. 
^We assume 8 hours per night throughout the analysis. 
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Figure 6: Sky coverage /sky of cosmic siiear survey as a function of i-band exposure time ioxp for 
different number of redshift bins; iV = 1 (red, solid), N = 2 (green, long-dashed) and = 3, 4 
(blue, short-dashed). Here, we fix the total observation time Tobs to 180 nights. The blue point on 
this line is located at Tobs = 180 nights where the survey parameters correspond to the values used 
in section 0. 
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Table 5: The number of redshift bin, N, redshift range for each redshift bin and ^^cxp for 
tomography. We use the i-band filter for the case = 1, and the r- and z-band filters for iV = 2. 
For the cases with = 3 and 4, we consider g-, r-, i- and z-band observations. 

5.2.1 Neutrino mass 

Let us first show the forecast results for the constraint on the total neutrino mass, fixing the 
total observation time to Tobs = 180 nights. In Figure ^, we plot the marginalized 1-a error, 
cr{Y^ rriy), as a function of the i-band exposure time texp- Here, the results are normalized 
by the value obtained from the canonical setup in section ^, i.e., (TM(^'mu) = 0.041 eV. 
Note again that all the power spectra are computed up to £max = 3000 to calculate the 
Fisher matrix. 

Basically, the resultant marginalized l-cr error is shown to be a monotonically increas- 
ing function of the exposure time. That is, the neutrino mass constraint is improved as 
decreasing the exposure time, and at texp ^ 3 min, the tightest constraint is obtained from 
the A'^ = 1 case. This indicates that for a limited observation time, a shallow and wide-field 
survey is preferable, and no tomographic technique is necessary to improve the neutrino 
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mass constraint. This result is contrasted with Hannestad et al. (2008) [|T5|, where they 
fixed the exposure time and sky coverage, and showed that the constraint on neutrino 
mass is improved with the number of redshift bins. However, the resultant survey setup for 
lensing tomography requires an increasingly long observation time, which seems in some 
sense impractical. 

A closer look at a shorter exposure time at t < 3 min reveals that the constraint with 
= 1 ceases to improve and turn next to gradually worsen. Eventually, the N = 2 case 
can give a better constraint on the neutrino mass. The basic reason for this behavior 
comes from the competition of the two effects. That is, as decreasing the exposure time, 
the survey area of cosmic shear experiment increases and the statistical error is reduced. 
However, as a trade off, the redshift distribution of galaxy samples becomes shallower and 
the number of samples per square arcminute tends to decrease, leading to the reduction of 
the signal-to-noise ratio. Hence, for N = 1 case, there appears an optimal exposure time, 
texp — 3 min, where the tightest constraint on the neutrino mass is obtained. 

Note that the optimal value of texp does also exist for A 7^ 1 cases, and these would be 
much shorter than that of the A = 1 case. In such a short exposure time with texp ^ 5 min, 
however, a robust and reliable shear measurement is practically difficult and challenging 
due to the stability of the seeing condition. Thus, we hereafter consider the exposure time 
of texp > 5 min, and look for the best constraint on the neutrino mass, characterized by 
the function of total observation time. 

Figure ^ shows the neutrino mass constraint plotted against the total observation time 
Tobs, which is determined by the optimal choice of the exposure time, so as to minimize the 
marginalized 1-a error a{^m^) for each observation time, under the condition tgxp ^ 5 
min. Again, the results are normalized by the fiducial value, crfid(J2^v)i obtained in the 
previous section. As anticipated from Figure 0, the optimal neutrino mass constraint is 
given as an decreasing function of Tobs, and roughly scales as T~^^'^. Under the condition 
^exp ^ 5 min, the A = 1 case can give the tightest constraint over the plotted range of Tobs- 
However, the resultant improvement factor, given by the inverse of mjy)/crfid(X] ?TT.,y), 
is rather small, and it somehow reaches at ~ 1.2 at Tobs = 250 nights. The dependence on 
the number of redshift bins is also very weak, and no significant improvement is expected 
from the optimization of the exposure time. 

5.2.2 Dark energy 

Next consider the constraints on dark energy properties, i.e., wq and Wa- In order to 
find the optimal exposure time for the dark energy equation-of-state parameters, we here 
introduce the figure-of-merit defined by @ 

FoM{wo,Wa) = , (5.4) 

VdetF-i 

where is the 2x2 sub-matrix whose components are taken from the inverse of Fisher 
matrix associated with the parameters wq and Wa- Since the determinant of F'~^ is pro- 
portional to the area of the two-dimensional error ellipse on wo-Wa plane, a larger value of 
the figure-of-merit implies a tighter constraint on dark energy properties. 
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Figure 7: The constraint on total neutrino mass, crQ2 a function of the exposure time texp 

for = 1 (red sohd), N = 2 (green long-dashed), N = 3 (blue short-dashed) and = 4 (magenta 
dotted). The total observation time Tobs is fixed to 180 nights. The resultant value of a(^TO^) is 
plotted dividing by the value obtained in section ^ (JRdiYl ^y)- 
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Figure 8: The expected neutrino mass constraint cr(X] '^i') ^ function of the total observation 
time Tobs, determined by the optimal choice of the exposure time. Here, the exposure time is 
optimized so as to minimize ai^mi,) for each observation time under the condition toxp > 5 min. 
The blue filled circle is the result for the fiducial survey setup; TV = 1 (red solid), TV = 2 (green 
long-dashed), TV = 3 (blue short-dashed) and TV = 4 (magenta dotted). 



In Figure |9|, fixing the total observation time to Tobs = 180 nights, we show the 
figure-of-merit as the function of exposure time, and the results are normalized by the 
value obtained from the canonical setup in section i.e., ¥oM^^{wo,Wa) = 168. The 
figure-of-merit is a decreasing function of exposure time, and among several choices of 
tomographic binning, the TV = 2 case can give a better performance for the figure-of-merit. 
While this result is consistent with Yamamoto et al. |57], the situation they considered is 
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Figure 9: The figure-of-merit, FoM{wo, Wa), as a function of exposure time, fixing total observa- 
tion 180 nights. The resuhs are normahzed by FoMgd(wo, ?«a)= 168, which is obtained from the 
canonical survey setup (see Sec.^). Line types are the same as in Figure 0. 

somewhat different from ours: they considered the cosmic shear experiment alone assuming 
the massless neutrinos. Also, the optimal number of filters is contrasted with the case of 
the neutrino mass constraint, i.e., = 1. Presumably, the differences between the optimal 
number of filters would reflect the fact that the dark energy properties sensitively affect 
the growth of the large-scale structure, while the effect of the massive neutrinos can give a 
free-streaming suppression on the small-scale structure growth, whose redshift dependence 
is extremely weak. In this sense, the tomographic technique of the cosmic shear experiment 
is especially helpful in constraining the time evolution of dark energy equation-of-state, and 
for our current setup with limited observation time, the choice of iV = 2 becomes optimal. 
Note that the choice of A = 2 still remains optimal even if we look for the maximum value 
of FoM including the neutrino mass. 

Next, in Figure IC, similarly to Figure 0, the figure-of-merit for various choices of A^ is 



plotted against the total observation time Tobs- Note here that varying texp, we looked for 
the maximum value of FoM for each observation time Tobs; under the condition iexp ^ 5 
min. Then, the best constraint is obtained from the A = 2 case, and the improvement 
of the FoM is roughly proportional to ~ Tobs- However, the resultant constraint is not so 
drastically changed compared to the one obtained from the canonical setup in section |^. 
Even at Tobs = 250 nights, the improvement is still below ~ 1.7. 

The results indicate that joint constraints on the neutrino mass and dark energy 
equation-of-state parameter is rather stable against the details of the survey setup of the 
cosmic shear experiment, and in this sense, we could say that the forecast results in section 
^ is robust. 



6. Summary 

In this paper, we explored a capability and a synergy of the two weak lensing experiments, 
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Figure 10: The figure-of-merit, FoM(u'o, Wa), as a function of total observation time Tobs, deter- 
mined by the optimal choice of the exposure time. Here, the exposure time is optimized so as to 
maximize FoM(i(;o, Wa) for each observation time under the condition ioxp > 5 min. The blue filled 
circle is the result for the fiducial survey setup. Line types are as the same in Figure 0. 



i.e., CMB lensing and cosmic shear, to simultaneously constrain the neutrino mass and 
dark energy properties. As representative lensing experiments, we consider the Planck 
and ACTPol for CMB lensing observations, and Subaru HSC survey for cosmic shear ex- 
periment. Including the primary CMB information as a prior cosmological information, 
the Fisher analysis showed that combining Planck, ACTPol, and HSC, the total mass of 
neutrinos with ~ 0.1 eV would be detected with a significance of ~ 2a level. As for the 
dark energy equation-of-state parameters, the combination of the Planck and HSC would 
provide a constraint on wq and Wa with the accuracy of 10% and 20% level, which roughly 
corresponds to the expected errors for the Stage-HI class experiment defined by the Dark 
Energy Task Force. In other words, including the neutrino mass in the parameter estima- 
tion would not drastically alter the Figure-of-Merit estimates of dark energy parameters in 
weak lensing measurements. The constraints will be further improved by a factor of ~ 2 if 
we add the ACTPol data. 

We have also studied the role of the cross correlation statistics obtained from the two 
lensing experiments. While the deflection-shear cross correlation gives a strong statistical 
correlation between CMB lensing and cosmic shear signals, which makes the improvement 
of the combined constraints rather moderate, the temperature-deflection cross correlation 
is found to be sensitive to the early-time evolution of dark energy, and it plays an important 
role to constrain Wa rather than t^o- 

Further, we have investigated the sensitivity to the choice of survey setup specifically 
focusing on the cosmic shear survey, and compare the results in the canonical setup (Sec. ^ 
with those for the optimal survey setup characterized by the exposure time and/or total 
observation time. Maximizing the efficiency of the survey setup while keeping the total 
observation time, the optimal number of redshift bins needed for the lensing tomography 
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becomes = 1 for neutrino mass and N = 2 for dark energy equation-of-state param- 
eters, which are contrasted with the canonical setup of = 3. However, the resultant 
improvement for the optimal survey setup is rather small, and the forecast results depend 
very weakly on the total observation time. Hence, for the setup with a limited observation 
time, forecast results for the canonical setup presented in section ^ seems robust against 
the details of the survey setup. 

Note, however, that this conclusion crucially depends on the assumption that the re- 
construction of the lensing deflection angle and cosmic shear measurement can be made 
perfectly without any serious systematics. Further, we assume that the theoretical tem- 
plate for lensing power spectra is well-understood, and can be described with the fitting 
formula [^. Also, the non-Gaussian covariance due to the nonlinear gravity has been 
neglected. These assumptions or treatment would be rather optimistic, and toward proper 
comparison with observations, the accurate template for lensing spectra should be devel- 
oped, and the effect of non-Gaussian covariance must be properly incorporated into the 
parameter estimation analysis [p4| , 51 1. In developing these issues, recently proposed an- 
alytical techniques to accurately predict matter power spectrum would be helpful (e.g., 
|11, pi, 16, 54, 41, |2^|), and a large set of N-body simulations is really demanding to study 



non-Gaussian covariance (e.g., pq|). 

Throughout the paper, we have focused on the homogeneous dark energy component 
which only affects the background dynamics of cosmic expansion. However, there may 
possibly exist a clustering component of dark energy (e.g. 1 20, |22[), which can alter the 
growth of structure on very large scales. As discussed by several authors (e.g.,||2^), CMB 
lensing and cosmic shear experiments are shown to be complementary probes for dark en- 
ergy clustering, and even including the effect of free-streaming of massive neutrinos, which 
appears on relatively small scales, they could still provide fruitful cosmological constraints 
on both neutrino and dark energy. In this respect, the scientific impacts from the two 
lensing experiments may become even more large. 

Finally, we note that combining the lensing data with other cosmological probes such 
as the baryon acoustic oscillations and type la supernovae can potentially break the degen- 
eracies between wq and Wa ||5^. Further, the clustering statistics of the galaxy distribution 
exhibit a free-streaming suppression by the massive neutrinos, and the joint analysis with 
lensing data would also break the degeneracy between ^ m^^ and il^h"^ ||4^, Although 
the control of the observational systematics such as galaxy bias and calibration of super- 
nova light curves become very much severe, the scientific impact on the neutrino mass and 
dark energy parameters would be significant, and a synergy between lensing and other 
cosmological probes should deserve further investigation. 
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A. Constraints on time varying dark energy equation-of-state parameter 

Section ^ reveals that the temperature-deflection cross correlation seems to play an im- 
portant role to constrain Wa, and it helps to break a degeneracy between Q\ and Wa- In 
this appendix, employing the principal component analysis (PCA) introduced by Ref. pSj , 
we investigate the sensitivity of the cross correlation statistics to the time evolution of 
dark energy in more general context, characterized by the time- varying equation-of-state 
parameter, w{z). Also, we explore the degeneracy between Q\ and w{z), and show that 
inclusion of the cross correlation helps to break the degeneracy even at high redshifts. 
To apply the PCA, let us we first discretize the equation-of-state parameter as 

N 

w{z) =^WiE{zi^i,Zi;z), (A.l) 

i=l 

where the function E{zi-i, Zi] z) is the step function defined by 1 at < z < Zi, and 
otherwise. Here, we set = 18 and choose the redshift interval as Zj — = 0.5. 
Treating the coefficients Wi as free parameters, we compute the statistical uncertainty of 
each coefficient from the Fisher matrix in the same way as described in section |^. Assuming 
the fiducial values of the coefficients Wi as —1 for i = I,-- - , A^, we consider the Fisher 
matrix for the lensing experiments combining Planck, ACTPol and HSC. Then, we obtain 
the N X N sub-matrix F for the coefficients Wi, marginalized over the other cosmological 
parameters'^. 

The resultant matrix F generically includes non-vanishing off-diagonal components, 
which represent statistical correlations between different Wi. That is, depending on the 
survey setup of the lensing experiments, there exist a strong degeneracy between different 
coefficients, and thereby the constraint on each coefficient Wi cannot be obtained inde- 
pendently. In this sense, a naive discretization of equation-of-state parameter ( [A.l| ) may 
not be a good description for time-varying equation-of-state parameter. Rather, a better 
characterization may be obtained by diagonalizing the matrix F: 

F = U^AU. (A.2) 

The quantity A is the diagonal matrix consisting of the eigenvalues of the matrix F, which 
we denote by Xj. The matrix U is the orthogonal matrix, whose element of i-th row and 
j-th column is given by the eigenvector ej{zi). Since the eigenvectors {ej{z)} can form an 

*To be precise, we first compute the inverse Fisher matrix, and construct the sub-matrix from it by 
extracting the components associated with the coefficients Wi. fnverting this sub-matrix again, we finaily 
get F. 
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orthonormal basis and can be regarded as function of z, we may expand the time-varying 
equation-of-state parameter as 

N 

w{z) = '^aiei{z) (A.3) 

i=l 

with the coefficients a, computed as ai = w{za)ei{za)- In this expansion, the coefficient 
ai is statisticaUy independent, and the associated error is inversely proportional to the 
eigenvalues Aj. In this sense, the functions ei{z) are the principal components (PCs), and 
the redshift dependence of the best-determined eigenvectors characterizes the sensitivity of 
the measurements to the time evolution of dark energy. 

Figure |Tl| shows three of the best-determined PCs for eigenvalues (left) and eigenvec- 
tors (right). Here, we separately plot the results in four cases, +Qd, +dd, +d, and +7 (see 
Sec. [4.2| and [4. 3D , among which the cosmic shear signal (+7) is shown to give the tightest 
constraint on the dark energy parameters. However, the redshift sensitivity of the cosmic 
shear is rather restricted to the lower redshift, and the eigenvectors become almost vanish- 
ing at z > 2. On the other hand, the eigenvectors associated with the CMB lensing signals 
(i.e., +d, +dd, +Qd) can have non- vanishing values even at higher redshifts. Although 
their eigenvalues are basically smaller than those of the cosmic shear, they still have some 
contributions to the constraints on the dark energy equation-of-state parameters, espe- 
cially at higher redshifts. An interesting point is that among several CMB lensing signals, 
temperature-deflection cross correlation seems to give an important contribution to the 
constraint on the early-time evolution of dark energy, since the first two largest eigenvalues 
of the temperature-deflection cross correlation are larger than those of the deflection auto 
correlation. 

To clarify the role of the temperature deflection cross correlation, we calculate the cor- 



relation coefficients between and Wi, r{il,\, Wi), the result of which is shown in Figure 12 



Note that the parameter Wi implies w{zi), and the resultant correlation coefficient is plotted 
against the redshift z. As increasing redshift, the absolute value of the correlation coeffi- 
cients for the temperature-deflection cross correlation becomes smaller than those for the 
deflection auto correlation. This indicates that the inclusion of temperature-deflection cross 
correlation helps to break a degeneracy between the density parameter and equation-of- 
state parameter, especially at higher redshifts. In this respect, the temperature-deflection 
cross correlation would be a more sensitive probe for the dark energy equation-of-state 
parameter than the deflection auto correlation. 
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